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In  Eq.  24,  7q  (u)  is  the  zeroth-order  Bessel  function  of  the  first  kind.  The  function  T  (^) 
is  the  complex  plane-wave  transmission  coefficient  and  is  related  to  the  reflection  coefficient 
of  Eq.  6  by 

r(^)  =  l  +  /?(^).  (29) 

This  transmission  coefficient  is  evaluated  at  the  complex  angle 
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1  -(-  K^-  2k  (sin0^.sin0^  +  cos0^cos0xos(j)^) 


This  angle  was  derived  by  assuming  that  the  local  tilt  of  the  surface  is  such  as  to  give 
Snell’s  law  transmission  in  the  scattered  direction.  This  is  done  in  analogy  with  the  in- water 
case,  where  assuming  that  the  local  tilt  of  the  surface  is  such  as  to  cause  specular  reflection  in 
the  scattered  direction  gives  results  closer  to  “exact”  numerical  calculations. 

Using  a  procedure  similar  to  that  of  Jackson  (1993),  one  can  approximate  the  result  in 
Eq.  24  as 
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The  term  10  (A^"^  +  lA^]"^)  /[A^P  in  Eq.  31  artificially  increases  the  Kirchhoff  contribution 
away  from  the  Snell’s  law  (refraction)  direction.  Since  Eq.  23  favors  the  smaller  of  the  Kirch¬ 
hoff  and  perturbation  approximations,  this  suppresses  the  Kirchhoff  contribution  except  near 
the  refraction  direction. 


3.1.2.  Perturbation  Theory  Cross  Section 

The  perturbation-theory  results  are  simplifications  of  those  presented  by  Moe  et  al. 
(1995).  In  particular,  the  sediment  below  the  rough  interface  is  assumed  to  be  nonlayered.  The 
resulting  cross  section  is 
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ABSTRACT 

A  model  is  developed  for  bistatic  scattering  into  ocean  sediments  by  an  acoustic  source 
located  in  the  water  and  a  receiver  buried  in  the  bottom.  The  model  is  closely  related  to  that 
developed  by  Jackson  (1993)  for  a  source  and  receiver  that  were  both  located  in  the  water  and 
parallels  that  effort  in  many  regards.  In  particular,  bottom  scattering  is  treated  as  being  due  to 
interface  roughness  and  volume  inhomogeneity.  Kirchhoff  theory  is  used  to  explain  interface 
scattering  near  the  refracted  (Snell’s  law)  direction,  and  perturbation  theory  (Moe  et  al.,  1995) 
is  used  for  all  other  directions.  Perturbation  theory  is  used  to  treat  the  scattering  due  to  volume 
inhomogeneities.  Model  results  are  presented  for  two  of  the  bottom  sites  examined  by  Jackson 
(1993),  and  possible  experimental  tests  of  the  model  are  discussed  briefly. 


TR  9505  iii 


UNIVERSITY  OF  WASHINGTON  •  APPLIED  PHYSICS  LABORATORY 


EXECUTIVE  SUMMARY 

This  report  describes  a  model  of  high-frequency  (10-100  kHz)  bistatic  bottom  scattering 
for  an  acoustic  source  located  in  the  water  column  and  a  receiver  buried  in  the  bottom  sedi¬ 
ment  (referred  to  here  as  in-sediment  scattering).  The  physics  incorporated  into  the  model  has 
much  in  common  with  previous  modeling  of  bistatic  scattering  done  at  the  Applied  Physics 
Laboratory,  University  of  Washington,  the  major  differences  being  that  in  the  previous  work 
both  the  source  and  the  receiver  were  located  in  the  water  column  (referred  to  here  as  the  in¬ 
water  case).  The  sediment  is  treated  as  a  lossy  fluid.  The  bistatic  scattering  is  hypothesized  to 
be  due  to  scattering  from  interface  roughness  and  from  volume  inhomogeneities  within  the 
sediment.  Interface  scattering  is  calculated  via  the  Kirchhoff  approximation  for  angles  near  the 
Snell’s  law  direction  and  via  perturbation  theory  for  all  other  directions.  Volume  perturbation 
theory  is  used  to  relate  bistatic  volume  scattering  to  three-dimensional  spectra  describing  fluc¬ 
tuations  in  sediment  density  and  compressibility. 

The  input  parameters  needed  (and  measured)  for  the  in-sediment  model  are  the  same  as 
those  used  for  the  in-water  case  with  the  addition  of  one  parameter  used  to  give  a  low- 
frequency  cutoff  to  the  spectra  of  the  volume  inhomogeneities.  This  low-wavenumber  cutoff 
is  required  to  prevent  volume-perturbation  theory  from  predicting  unrealistically  large  scatter¬ 
ing  near  the  forward-scattering  direction.  The  other  parameters  are  the  sediment/water  sound- 
speed  ratio,  the  sediment/water  density  ratio,  two  parameters  that  define  the  two-dimensional 
interface  relief  spectrum,  and  three  parameters  that  characterize  the  three-dimensional  spectra 
describing  fluctuations  in  sediment  density  and  compressibility. 

One  of  the  attributes  of  the  model  is  that  it  predicts  penetration  of  acoustic  energy  into  the 
sediment  at  “post-critical”  angles,  i.e.,  grazing  angles  less  than  the  critical  angle  predicted  for 
a  sediment  with  a  sound  speed  greater  than  that  in  water.  There  has  been  evidence  of  post- 
critical  penetration  for  several  years  (Chotiros,  1995),  but  the  mechanism  involved  is  still  a 
subject  of  research.  In  this  model,  diffraction,  as  calculated  via  surface-interface-perturbation 
theory,  is  responsible  for  most  of  this  post-critical  penetration.  In  fact,  the  perturbation  model 
adopted  here  has  been  used  to  examine  the  experimental  evidence  of  post-critical  penetration 
(Moe  et  al.,  1995).  The  conclusion  at  this  point  is  that  this  mechanism  predicts  results  consis¬ 
tent  with  those  experiments. 
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1.  INTRODUCTION 

The  model  presented  here  borrows  heavily  from  the  bistatic  model  of  Jackson  (1993). 
Both  the  physical  mechanisms  incorporated  into  the  model  and  the  environmental  parameters 
needed  as  model  inputs  are  the  same.  The  difference  is  that  here  the  goal  is  to  determine  the 
scattering  into  the  sediment  (referred  to  as  in-sediment  bistatic  scattering)  whereas  Jackson 
was  interested  in  determining  the  scattering  into  the  water  column  (referred  to  as  in-water  bi¬ 
static  scattering). 

The  variables  and  form  of  the  equations  have  been  chosen  so  as  to  closely  parallel  the  pre¬ 
vious,  in-water  case.  However,  the  detailed  expressions  for  any  particular  variable  may  be  dif¬ 
ferent.  Section  2  contains  some  basic  definitions  and  defines  the  input  parameters.  Section  3 
gives  the  expressions  for  both  the  roughness  and  the  volume-scattering  portions  of  the  model. 
Section  4  presents  sample  results  for  two  sites  (Quinault  and  Arafura  Sea)  that  were  also  ex¬ 
amined  by  Jackson  (1993)  and  briefly  discusses  an  experimental  test  of  the  model  that  could 
be  carried  out.  Section  5  is  a  summary. 

2.  DEFINITIONS 

The  bistatic  scattering  strength  will  be  written  in  the  form 

h  't’s'  =  'Olog  (9^,  (1,^,  9,,)  +  (9,,  9,)  ]  ,  (1) 

where  log  is  base  10  and  (0^,  (f)^,  6.)  and  (0^,  (j)^,  0,)  are,  respectively,  the  rough¬ 
ness  and  volume  contributions  to  the  scattering  cross  section  per  unit  area.  The  definition  of 
scattering  cross  section  implicit  in  the  results  given  here  is  slightly  different  than  that  conven¬ 
tionally  used.  This  definition  is  given  in  the  appendix. 

2.1  Geometry 

The  angles  0^ ,  (j)^ ,  and  0-  in  Eq.  1  are  defined  in  Figure  1;  0^  is  the  incident  grazing  an¬ 
gle,  0^  is  the  scattered  grazing  angle,  and  is  the  bistatic  angle  defined  as  the  difference  in 
azimuth  between  the  incident  and  scattered  angles.  In  general,  one  needs  four  angles,  two 
grazing  angles  and  two  azimuthal  angles,  but  only  the  azimuthal  difference  is  needed  here  be¬ 
cause  bottom  statistics  are  assumed  to  be  transversely  isotropic. 
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Figure  1.  Bistatic  scattering  geometry.  The  incident  grazing  angle  is  denoted  0.,  the 
scattered  grazing  angle  is  denoted  0^,  and  the  bistatic  angle,  defined  as  the 
difference  in  azimuth  between  the  incident  and  scattered  directions,  is  denoted  . 


The  vectors  k.  and  are  also  shown  in  Figure  1.  These  vectors  are  the  incident  and  scat¬ 
tered  wave  vectors,  respectively.  From  Figure  1,  we  have 

k-  =  k{xcosQ-~zsmd-)  (2) 

and 

k^=  k^[x  {cosQ^cos(^^)  -i- j  (cos0^sin(t)^)  -zsin0J  ,  (3) 

where  k  is  the  acoustic  wavenumber  in  the  water  and  is  the  (complex)  acoustic  wavenum¬ 
ber  in  the  sediment. 
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2.2  Model  Inputs 

The  input  parameters  for  the  bistatic  bottom-scattering  model  are  given  in  Table  1.  Pro¬ 
cedures  have  been  developed  to  provide  estimates  of  these  parameters  when  direct  physical 
measurements  are  lacking  (Mourad  and  Jackson,  1989).  The  bottom  material  is  treated  as  a 
fluid  with  no  layering  or  gradients  (apart  from  the  random  fluctuations  responsible  for  volume 
scattering).  The  density  ratio,  sound-speed  ratio,  and  loss  parameter  characterize  the  nonran¬ 
dom  acoustic  properties  of  the  sediment.  It  is  convenient  to  define  the  following  complex 
quantities  in  terms  of  these  parameters  (Mourad  and  Jackson,  1989). 


K 


1  +  /5 


(4) 


P(0)  =  Vk^- cos^0  , 


(5) 


Table  1.  Model  input  parameters.  Also  needed  are  values  for  frequency,  water  sound  speed, 
and  water-mass  density. 


Symbol 

Definition 

Short  Name 

P 

Ratio  of  sediment  mass  density  to 
water  mass  density 

Density 

ratio 

V 

Ratio  of  sediment  sound  speed  to 
water  sound  speed 

Sound-speed 

ratio 

5 

Ratio  of  imaginary  wavenumber  to 
real  wavenumber  for  the  sediment 

Loss 

parameter 

Ys 

Exponent  of  sediment  inhomogeneity 
spectrum 

Inhomogeneity 

exponent 

^3 

Strength  of  sediment  inhomogeneity  spectrum 
(cm^)  at  wavenumber  {2k} /X  =  1  cm"* 

Inhomogeneity 

strength 

P 

Ratio  of  compressibUity  to  density 
fluctuations  in  the  sediment 

Fluctuation 

ratio 

Y2 

Exponent  of  the  bottom  relief  spectrum 

Spectral 

exponent 

^2  . 

Strength  of  bottom  relief  spectrum  (cm'*) 
at  wavenumber  {2k} /X  =  1  cm"* 

Spectral 

strength 

a 

Low-wavenumber  cutoff  parameter  for  sediment 
inhomogeneity  spectrum  (cm) 

Cutoff 

parameter 
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and 


^  (psine)/[P(e)]  -  1 
(psine)/[P(e)] +r 


(6) 


In  Eq.  5,  the  quantity  inside  the  square  root  is  complex,  and  the  square  root  is  defined  as 
the  complex  number  whose  magnitude  is  the  positive  square  root  of  the  original  magnitude 
and  whose  angle  is  one-half  the  original  angle.  The  quantity  P  (6),  defined  by  Eq.  5,  is  the 
normalized  z  component  of  the  wave  vector  in  the  sediment  for  an  incoming  plane  wave  hav¬ 
ing  grazing  angle  0  in  the  water.  This  definition  assumes  that  the  bottom  is  perfectly  flat  and 
homogeneous,  although  the  scattering  models  to  be  used  later  do  not.  The  normalization  is 
simply  division  by  k,  the  acoustic  wavenumber  in  the  water.  Finally,  P  (6)  is  the  flat-interface 
Rayleigh  reflection  coefficient.  The  complex  quantities  P  (0)  and  P  (0)  will  appear  in  later 
expressions  for  the  roughness  and  sediment  volume-scattering  cross  sections.  The  quantity  k 
is  the  ratio  of  the  acoustic  wavenumber  in  the  sediment  (a  complex  number)  to  the  acoustic 
wavenumber  in  the  water  (a  real  number).  The  loss  parameter,  5 ,  can  be  obtained  in  terms  of 
the  more  commonly  defined  sediment  sound  absorption  coefficient,  ,  as  follows: 


a^vclnlO 

40k/ 


(7) 


Care  must  be  exercised  in  the  choice  of  units  here.  If  is  in  units  of  decibels  per  meter, 
as  is  usually  the  case,  and  the  water  sound  speed,  c,  is  in  meters  per  second,  then  the  acoustic 
frequency,/,  must  be  in  hertz. 


The  density  ratio,  sound-speed  ratio,  and  loss  parameters  are  identical  to  those  used  by 
Mourad  and  Jackson  (1989)  in  their  backscatter  model,  as  are  the  two  parameters  used  to  de¬ 
scribe  random  bottom  relief  (also  referred  to  as  bottom  roughness).  Bottom  relief  is  assumed 
to  be  a  Gaussian  random  process  having  the  following  two-dimensional  spectrum: 


^2 

W(P)  =  - 1—.  (8) 

In  this  expression,  Wj  and  ^2  ^^e  model  parameters  for  the  strength  and  exponent  of  the 
relief  spectrum,  and  is  simply  a  reference  length  equal  to  1  cm.  The  subscript  “2”  denotes 
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the  fact  that  bottom  relief  is  a  two-dimensional  random  process.  Since  the  relief  spectrum  de¬ 
pends  only  on  the  magnitude  of  the  two-dimensional  wave  vector,  {K  =  \K\ ),  bottom  relief 
is  assumed  to  be  isotropic.  We  will  use  upper-case  letters  (e.g.,  K)  to  denote  transverse  vec¬ 
tors.  The  spectrum  is  normalized  so  that 


2k  j  W  (K)  KdK  =  (K^) .  (9) 

K 

C 

In  this  expression,  h'^  (K^)  is  the  mean- square  roughness  due  to  those  Fourier  compo¬ 
nents  of  bottom  relief  having  a  wavenumber  greater  than  .  The  wavenumber  can  be 
viewed  as  a  cutoff  on  the  spectrum.  It  is  an  independent  variable  employed  to  yield  finite  val¬ 
ues  for  mean-square  roughness.  The  spectral  exponent  is  constrained  as  follows: 

2<Y2<4.  (10) 

The  relief  statistics  can  be  alternately  described  by  the  “structure  function”  (Yaglom, 
1962),  itself  a  power  law  if  the  spectrum  is  a  power  law: 

DiR)  =  ([/(/Jo  +  i?)-/(/?o)]2)  =  c2/?2a,  (H) 

where/  (R)  is  the  z  coordinate  of  the  interface  at  horizontal  position  R  =  (x,  y).  The  func¬ 
tion  f  (R)  is  a  zero-mean,  two-dimensional  Gaussian  random  process.  The  structure  function 
parameters  are  related  to  the  spectrum  parameters  as  follows: 


a  =  y-l  (12) 

27tw.,r(2-a)2-2a 

- f - ,  (13) 

;zj2a(l  -a)r(l  -Fa) 

where  T  is  the  gamma  function. 

The  four  new  parameters  (jy  Wy  M-’  1^®  in-sediment  (and  in-water)  bistatic 

model  but  not  in  earlier  monostatic  models  describe  three-dimensional  sediment  inhomogene¬ 
ity.  These  replace  a  single  volume- scattering  parameter  appearing  in  the  older  monostatic 


and 


C2 
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models.  Sediment  inhomogeneities  are  described  by  three-dimensional  spectra  for  compress¬ 
ibility  and  density  fluctuations,  normalized  by  dividing  by  their  mean  values.  We  denote  these 
normalized  compressibility  and  density  fluctuations  by  (r)  and  Yp  (r)  ,  respectively.  We 
assume  that  the  fluctuations  in  sediment  compressibility  and  density  are  tied  together  by  a  sim¬ 
ple  proportionality  relation, 


=  liyp(r)  .  (14) 

In  saturated  sediments,  we  expect  the  model  parameter  p  to  have  values  in  the  neighbor¬ 
hood  of  -1  because  density  and  compressibility  tend  to  be  anticorrelated.  For  example,  if  a 
loosening  of  the  sediment  produces  a  10%  decrease  in  density,  then  there  will  be  a  correspond¬ 
ing  increase  of  about  10%  in  compressibility.  These  two  changes  must  be  of  about  the  same 
magnitude,  because  it  is  known  that  a  10%  change  in  density  is  likely  to  produce  a  much  small¬ 
er  percentage  change  in  sediment  sound  speed,  .  Since 


(15) 


it  follows  that  density  and  compressibility  changes  must  be  anticorrelated  in  order  to  produce 
a  small  change  in  sound  speed.  (See  Hamilton  and  Bachman  (1982)  and  Richardson  and 
Briggs  (1993)  for  data  on  the  relationship  between  density  and  sound  speed.) 


We  assume  that  the  fluctuations  in  density  and  compressibihty  are  spatially  stationary,  so 
that  power  spectra  can  be  defined  as  Fourier  transforms  of  the  corresponding  fluctuation  co- 
variances. 
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We  assume  that  the  spectrum  for  density  fluctuations  obeys  the  following  power  law: 

{1.0-exp[-(t^a)2/2]}  .  (19) 

Equation  14  then  forces  the  following  expressions  for  the  compressibility  fluctuation  spec¬ 
trum,  (k^),  and  the  cross  spectrum, 

2 

=  - ^{l-0-exp[-(/:^a)V2]}  (20) 

and 


=  - ^-{1.0-exp[-(l;  p)V2]}\  (21) 

For  the  sake  of  simplicity,  three-dimensional  isotropy  is  assumed  for  the  sediment  inho¬ 
mogeneity.  This  is  reflected  in  Eqs.  16-21  in  that  these  spectra  depend  only  upon  the  magni¬ 
tude,  A-y ,  of  the  three-dimensional  wave  vector,  Layered  sediments,  by  definition,  do  not 
possess  three-dimensional  isotropy.  It  is  an  open  question  whether  a  more  complicated,  aniso¬ 
tropic  sediment  model  is  needed  at  high  frequencies. 


The  term  in  braces  in  Eqs.  19-21  introduces  a  low  cutoff  into  the  spectra.  This  pre¬ 
vents  an  unphysical  divergence  near  forward  scattering.  In  the  results  discussed  in  this  report, 
a  is  set  to  50  cm.  The  filtered  spectra  of  Eqs.  19-21  then  agree  with  the  unfiltered  ones 
(Eqs.  19-21  without  the  term  in  braces)  down  to  0(1  m).  Furthermore,  integration  of  Eq.  19 
(for  the  case  =  3)  gives 


471^3  j£/[-(A2fl2)] 
2 


+  ln(A^) 


(22) 
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for  the  density  variance,  where  is  a  high-wavenumber  cutoff  such  that  the  smallest  vol¬ 
ume  inhomogeneities  are  0(1  mm)  and  Ei  is  the  exponential  integral  function.  When  this 
result  is  evaluated  using  a  =  50  cm  and  the  values  for  the  sediments  in  Section  4,  the  stan¬ 
dard  deviation  is  less  than  0.18,  which  is  consistent  with  available  core  density  data  (Jackson, 
1993). 


3.  BISTATIC  SCATTERING  CROSS  SECTIONS 
3.1  Surface  Roughness  Scattering 

In  this  section,  expressions  are  presented  for  the  bistatic  cross  sections  in  the  Kirchhoff 
and  perturbation  approximations.  The  net  cross  section,  (0^,  (j)^,  0.),  appearing  in  Eq.  1  is 
formed  by  smooth  interpolation  between  the  Kirchhoff  cross  section  near  the  specular  direc¬ 
tion  and  the  perturbation-theory  cross  section  elsewhere.  The  interpolation  between  these  two 
approximations  is  defined  so  that  the  smaller  cross  section  takes  precedence.  This  procedure 
is  based  on  the  fact  that,  for  power-law  spectra,  the  perturbation  approximation  overpredicts 
scattering  near  the  specular  direction  owing  to  the  singularity  in  the  relief  spectrum  at  zero 
wavenumber.  In  contrast,  the  Kirchhoff  approximation  tends  to  overpredict  in  other  direc¬ 
tions.  The  interpolation  scheme  used  here  is 

with  T|  =  -2 .  The  expressions  for  the  Kirchhoff  (a^^p  and  perturbation-theory  (C7^^)  cross 
sections  are  defined  next. 


3.1.1.  Kirchhojf  Cross  Section 

In  analogy  with  the  in-water  monostatic  expression  used  by  Jackson  and  Briggs  (1992) 
and  the  bistatic  expression  used  by  Jackson  (1993),  the  Kirchhoff  bistatic  cross  section  is 


|r(^)|2 

Ia^I  T 

87t 

A  A, 

\  z  tj 

je  ^^^^J^{u)udu  , 
0 


(24) 


where 
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(^|A  1)2 

<7  =  — (25) 

2 

=  {(-K  +  Ksin0-sin6^  +  Kcos0^cos6-cos(t)^) 

(2—T')  ’  (26) 

+  pA-y-^  (1  -  Ksin0.sin0_j,-Kcos0^cos0jCOs0^)  } 

=  (Ksin0^- sin0p  ,  (27) 

A  =  fcos20.--cos0.cos0  cos(j)  +4;cos^0  (28) 

1  i  V  '  ^ 

InEq.  24,  /q  («)  is  the  zeroth-order  Bessel  function  of  the  first  kind.  The  function  T  (^) 
is  the  complex  plane-wave  transmission  coefficient  and  is  related  to  the  reflection  coefficient 
of  Eq.  6  by 


Ta)  =  i-F/?(^). 


This  transmission  coefficient  is  evaluated  at  the  complex  angle 


f 

^  =  acos 

V 


k2  [  l-(sin0.sin0  -F  COS0  cos0,cos(})  )  2]  1  6-5^ 

_ *  ^ _ £ _ i _ £ _ > 

1  +  K2-2K(sin0j.sin0_^-F  cos0^cos0.cos(j)p  J  ^ 


(29) 


(30) 


This  angle  was  derived  by  assuming  that  the  local  tilt  of  the  surface  is  such  as  to  give 
Snell’s  law  transmission  in  the  scattered  direction.  This  is  done  in  analogy  with  the  in-water 
case,  where  assuming  that  the  local  tilt  of  the  surface  is  such  as  to  cause  specular  reflection  in 
the  scattered  direction  gives  results  closer  to  “exact”  numerical  calculations. 


Using  a  procedure  similar  to  that  of  Jackson  (1993),  one  can  approximate  the  result  in 
Eq.  24  as 

^  *  Ai  A  |4x  \  1 


,) 


1  + 10 


IaT 


(A,‘'“  +  ai72|Aj‘') 


1  -(-  a 
2a 


(31) 
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where 


(32) 


The  term  10  (A^  +  |A3|4)  /|a2|2  in  Eq.  31  artificially  increases  the  Kirchhoff  contribution 
away  from  the  Snell’s  law  (refraction)  direction.  Since  Eq.  23  favors  the  smaller  of  the  Kirch¬ 
hoff  and  perturbation  approximations,  this  suppresses  the  Kirchhoff  contribution  except  near 
the  refraction  direction. 

3.1.2.  Perturbation  Theory  Cross  Section 

The  perturbation-theory  results  are  simplifications  of  those  presented  by  Moe  et  al. 
(1995).  In  particular,  the  sediment  below  the  rough  interface  is  assumed  to  be  nonlayered.  The 
resulting  cross  section  is 

=  Aik,Q^)]^a^[l+R{d.)]  +foj[l-i?(ep]|V(A^p,  (35) 

where 


A(k,  ep 


(36) 
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1  cos^G„ 

'l - 

1  cos^G  /  cos^G 

Al 

1  \cos6  cos6  cos(t)„  „2 

t  S  S  ^  iv 

fli  =  I  -  - 1  - +  1  — , 

J  V  p 


= 


cos^G. 


fl- 


v2  y 


sinGpi  -p)  , 


(37) 


(38) 


(39) 


=  ^(A^  + 0.0001)  .  (40) 

The  expression  for  W  is  given  in  Eq.  8  and  that  for  A^  in  Eq.  28.  The  definition  in  Eq.  40 
is  used  to  avoid  the  singularity  in  the  relief  spectrum  at  zero  wavenumber.  In  practice,  the  di¬ 
vergence  is  not  a  problem  since  in  this  regime  the  Kirchhoff  result  is  used.  The  term  0.0001  in 
Eq.  40  simply  keeps  computer  codes  from  self-destructing. 


3.2  Volume  Scattering 

As  in  the  in-water  case,  an  effective  interface-scattering  cross  section  is  used  that  is  de¬ 
rived  from  an  integration  of  the  volume  scattering  in  depth.  As  a  note  of  caution,  this  proce¬ 
dure  may  lead  to  an  overprediction  of  volume  scattering  near  normal  incidence.  In  practice, 
however,  this  overprediction  will  probably  not  be  visible,  since,  near  normal  incidence,  the 
Kirchhoff  interface  scattering  is  the  dominant  process  in  getting  energy  to  the  in-sediment  re¬ 
ceiver.  (See  Section  4.) 

Using  conventional  volume  scattering  perturbation  theory,  Jackson  (1993)  arrived  at  an 
expression  for  the  volume-scattering  cross  section.  The  same  result  applies  here. 
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4.  MODEL  ILLUSTRATION  AND  EXPERIMENTAL  TESTS 

We  illustrate  the  model  by  employing  parameters  appropriate  to  two  sites,  Quinault  and 
Arafura  Sea.  The  Quinault  site  has  a  sandy  bottom  and  is  off  the  coast  of  Washington  State. 
The  Arafura  Sea  site  has  a  silty  bottom  and  is  north  of  Australia.  Table  2  gives  the  model  pa¬ 
rameters  for  both  sites. 


Table  2.  Model  input  parameters  for  Quinault  and  Arafura  Sea  sites. 


Site 

Sound- 

Speed 

Ratio 

(V) 

Density 

Ratio 

(P) 

Loss 

Param. 

(6) 

Inhomo. 

Expon. 

(Yy) 

Inhomo. 

Strength 

(Wo), 

3 

cm 

Fluct. 

Ratio 

(|I) 

Cutoff 

Parameter 

(a), 

cm 

Relief 

Expon. 

(Y2) 

Relief 

Strength 

(W,), 

cm 

0.01150 

3.0 

0.000127 

-1.0 

50 

3.67 

0.00422 

0.9880 

1.390 

0.00533 

3.0 

-1.0 

50 

3.18 

0.00318 

For  the  these  sites,  all  parameters  except  those  for  volume  inhomogeneity  are  taken  from 
Jackson  and  Briggs  (1992).  The  inhomogeneity  exponent,  ,  was  set  to  3.  As  noted  by  Jack- 
son  and  Briggs,  this  yields  a  frequency-independent  contribution  of  volume  scattering  to  the 
effective  interface  scattering  strength,  in  agreement  with  available  data.  The  parameter  |i , 
which  relates  density  and  compressibiUty  fluctuations,  is  taken  to  be  -1.  The  inhomogeneity 
strength  parameter,  Wg ,  was  determined  by  requiring  that  the  normalized  volume-scattering 
parameter  (Eq.  46  of  Jackson,  1993)  have  the  value  used  by  Jackson  and  Briggs.  As  mentioned 
earher,  the  interpolation  parameter,  Tj ,  was  set  to  -2. 

4.1  Scattering  Strength  vs  Angle 

Figures  2-7  show  predictions  of  the  total  bistatic  scattering  strength  and  the  scattering 
strength  associated  with  each  component  of  the  model  as  a  function  of  scattered  grazing  angle, 
with  the  incident  grazing  angle  and  bistatic  angle  held  fixed.  Figures  2-4  are  for  the  Quinault 
site  and  Figures  5-7  are  for  the  Arafura  site.  The  goal  of  this  section  is  both  to  give  results  that 
can  be  used  to  check  code  implementation  and  to  comment  on  some  of  the  characteristics  seen 
in  the  figures. 
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Figure  2.  Model  predictions  of  bistatic  scattering  strength  at  40  kHz  for  Quinault  site 
as  a  function  of  scattered  grazing  angle,  with  the  incident  grazing  angle 
fixed  at  90°  and  the  bistatic  angle  fixed  at  180°.  Separate  components  and 
the  interpolated  total  scattering  are  shown. 


Scattered  Grazing  Angle  (deg) 


-  Total 

. .  Kirchhoff 

- Surface  Pert. 

- Volume  Pert. 


Figure  3.  Model  predictions  of  bistatic  scattering  strength  at  40  kHz  for  Quinault  site 
as  a  function  of  scattered  grazing  angle,  with  incident  grazing  angle  fixed 
at  45°  and  bistatic  angle  fixed  at  180°.  Separate  components  and  the 
interpolated  total  scattering  are  shown. 
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Figure  4.  Model  predictions  of  bistatic  scattering  strength  at  40  kHz  for  Quinault  site 
as  a  function  of  scattered  grazing  angle,  with  incident  grazing  angle  fixed 
at  10°  and  bistatic  angle  fixed  at  180°  Separate  components  and  the 
interpolated  total  scattering  are  shown. 
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- Surface  Pert. 
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Figure  5.  Model  predictions  of  bistatic  scattering  strength  at  40  kHz  for  Arafura  site 
as  a  function  of  scattered  grazing  angle,  with  incident  grazing  angle  fixed 
at  90°  and  bistatic  angle  fixed  at  180°.  Separate  components  and  the 
interpolated  total  scattering  are  shown. 
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Figure  6.  Model  predictions  of  bistatic  scattering  strength  at  40  kHz  for  Arafura  site 
as  a  function  of  scattered  grazing  angle,  with  incident  grazing  angle  fixed 
at  45°  and  bistatic  angle  fixed  at  180°.  Separate  components  and  the 
interpolated  total  scattering  are  shown. 
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Figure  7.  Model  predictions  of  bistatic  scattering  strength  at  40  kHz  for  Arafura  site 
as  a  function  of  scattered  grazing  angle,  with  incident  grazing  angle  fixed 
at  10°  and  bistatic  angle  fixed  at  180°.  Separate  components  and  the 
interpolated  total  scattering  are  shown. 
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When  examining  the  figures,  it  is  important  to  remember  that  the  interpolation  scheme  of 
Eq.  23  results  in  the  smaller  of  either  the  Kirchhoff  or  the  perturbation  results  being  used  in 
calculating  the  total  scattering  strength.  For  the  Quinault  site,  the  sound-speed  ratio  in  Table  2 
implies  the  existence  of  a  critical  angle  of  about  26®.  Therefore  Figures  2  and  3  are  for  incident 
angles  where  a  refracted  path  would  be  expected  for  a  flat  surface.  The  peak  in  the  scattering 
strength  is  at  the  scattered  angle  appropriate  for  a  refracted  path,  and  the  Kirchhoff  approxi¬ 
mation  is  being  used  to  calculate  the  scattering  strength  near  this  (Snell’s  law)  direction.  Fur¬ 
ther  away  from  this  direction,  the  surface-perturbation  contribution  first  comes  into  play  and 
then  the  volume-scattering. 

The  volume-scattering  contribution  has  been  approximated  in  terms  of  an  equivalent 
surface-scattering  cross  section  (given  by  Eq.  45);  i.e.,  the  denominator  in  Eq.  45  is  a  measure 
of  the  relative  contribution  of  volume  scatterers  as  a  function  of  depth.  In  Figure  2,  where  the 
incident  grazing  angle  is  90®,  the  denominator  becomes  very  small  near  scattered  grazing  an¬ 
gles  of  90®.  Physically,  this  is  because  the  increased  attenuation  of  the  incident  wave,  in  pro¬ 
pagation  from  the  surface  to  the  scattering  volume,  is  balanced  by  a  decrease  in  attenuation  in 
propagation  from  the  scattering  volume  to  the  receiver.  As  a  result,  the  present  approximation 
almost  certainly  overpredicts  the  volume-scattering  strength  for  these  near-normal-incidence 
geometries.  However,  since  the  surface  roughness  (the  Kirchhoff  contribution)  is  the  dominat¬ 
ing  effect  for  scattered  grazing  angles  near  90®,  this  overprediction  would  seem  to  be  of  little 
practical  significance.  Nevertheless,  caution  is  warranted  when  using  the  present  model  near 
normal  incidence. 

Figure  4  is  for  an  incident  angle  well  past  critical,  so  a  refracted  path  into  the  sediment 
does  not  exist  and  the  Kirchhoff  approximation  never  plays  a  role.  The  volume-scattering  con¬ 
tribution  (which  results  from  an  evanescent  wave  in  the  sediment  scattering  from  volume  in¬ 
homogeneities)  is  smaller  than  the  surface-perturbation  contribution  for  all  but  a  small  range 
of  scattered  angles  near  backscattering.  The  results  given  in  Section  4.2  will  show  that  the 
scattered  grazing  angles  near  90®  play  an  important  role  in  the  signal  that  a  buried  receiver  will 
see.  Near  90,®  surface  perturbation  theory,  i.e.,  diffraction  due  to  surface  roughness,  domi¬ 
nates  the  scattering  process. 

For  the  Arafura  site,  the  sound  speed  in  the  sediment  is  less  than  that  in  the  water,  so  no 
critical  angle  exists.  Figures  5-7  show  that  the  scattering  strength  peaks  in  direction  of  the 
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refracted  path.  Furthermore,  the  surface  roughness  leads  to  a  much  narrower  peak  in  this  di¬ 
rection.  This  will  be  seen  to  have  ramifications  on  the  figures  shown  in  the  following  section. 
Finally,  volume  scattering  plays  the  major  role  in  the  scattering  strength  away  from  this  nar¬ 
row  peak.  The  coloring  scheme  used  in  the  figures  masks  the  blue  of  volume  scattering  by  the 
red  of  total  scattering  when  the  two  are  nearly  equal. 

4.2  Scattering  Strength  vs  Position 

Typically,  the  scattering  strength  will  be  used  to  calculate  the  intensity  received  at  a  point 
in  the  sediment  from  a  source  in  the  water  column.  This  involves  an  integral  over  the  water/ 
sediment  interface  in  which  the  integrand  includes  the  scattering  cross  section  at  each  interface 
position,  spherical  spreading  from  the  source  to  the  surface  point  then  from  the  surface  point 
to  the  receiver,  and  attenuation  within  the  sediment.  Figures  8-21  present  plots  of  (1)  the  scat¬ 
tering  strength  as  a  function  of  position,  and  (2)  the  integrand  in  the  calculation  of  received 
intensity  for  three  source/receiver  geometries.  The  scattering-strength  plots  indicate  the  rela¬ 
tive  levels  of  intensity  that  traverse  the  interface,  whereas  the  integrand  plots  show  how  these 
scattering  strengths  contribute  to  the  received  intensity. 

In  all  the  geometries  examined,  the  horizontal  separation  between  the  receiver  and  the 
point  source  is  20  m,  and  the  depth  of  the  receiver  below  the  water/sediment  interface  is  1  m. 
The  heights  of  the  source  above  the  interface  are  5, 10,  and  20  m.  The  source  and  receiver  po¬ 
sitions  are  given  on  each  figure  in  x,y,z  format.  Figures  8-13  are  for  the  Quinault  site,  and 
Figures  14-21  for  the  Arafura  site. 

Figures  8,10,  and  12  show  the  scattered  strength  as  a  function  of  position  on  the  interface 
for  source  heights  of  20,  10,  and  5  m,  respectively.  As  the  source  is  lowered,  the  elliptical  re¬ 
gion  where  the  highest  intensity  is  scattered  into  the  sediment  elongates,  and  the  peak  intensity 
moves  toward  the  source.  The  position  of  the  peak  corresponds  to  the  refracted  path  for  a  flat 
interface.  There  is  also  a  secondary  maximum  in  the  scattering  strength  immediately  below 
the  source.  This  is  due  to  volume  scattering.  Contrast  these  figures  with  Figures  9, 1 1,  and  13, 
which  show  the  relative  contributions  to  the  intensity  seen  at  the  receiver  for  source  heights  of 
20,  10,  and  5  m.  The  most  significant  contributions  come  from  the  region  immediately  above 
the  receiver.  This  is  because  sound  from  that  region  has  the  shortest  path  to  traverse  within  the 
sediment  (and  therefore  the  least  attenuation)  to  arrive  at  the  receiver. 
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Figure  8.  Scattering  strength  at  40  kHz  as  a  function  of  position  for  source  and 
receiver  positions  given  in  meters  at  left  (x,y,z),  Quinault  site.  For  this 
geometry,  the  incident  grazing  angle  immediately  above  the  receiver,  i.e., 
at  (20,0,0),  is  well  above  the  critical  angle  of  26°. 
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Figure  9.  Scattering  strength  of  Figure  8  minus  (see  Eq.  7)  times  distance 
traveled  from  each  surface  position  to  the  receiver.  The  result  indicates  the 
major  region  of  the  surface  contributing  to  the  received  energy.  All  values 
below  -90  dB  have  been  set  equal  to  -90  dB.  Note  that  the  majority  of  the 
plot  is  at  the  red  level  corresponding  to  -90  dB. 
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Figure  10.  Scattering  strength  at  40  kHz  as  a  function  of  position  for  source  and 
receiver  positions  given  in  meters  at  left  (x,y,z),  Quinault  site.  For  this 
geometry,  the  incident  grazing  angle  of  26.5°  immediately  above  the 
receiver,  i.e.,  at  (20,0,0),  is  approximately  equal  to  the  critical  angle  of 
26°. 
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Figure  11.  Scattering  strength  of  Figure  10  minus  (see  Eq.  7)  times  distance 
traveled  from  each  surface  position  to  the  receiver.  The  result  indicates 
the  major  region  of  the  surface  contributing  to  the  received  energy.  All 
values  below  -90  dB  have  been  set  equal  to  -90  dB. 
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Figure  12.  Scattering  strength  at  40  kHz  as  a  function  of  position  for  source  and 
receiver  positions  given  in  meters  at  left  (x,y,z),  Quinault  site.  For  this 
geometry,  the  incident  grazing  angle  of  14°  immediately  above  the 
receiver,  i.e.,  at  (20,0,0),  is  well  below  the  critical  angle  of  26°. 
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Figure  13.  Scattering  strength  of  Figure  12  minus  (see  Eq.  7)  times  distance 
traveled  from  each  surface  position  to  the  receiver.  The  result  indicates 
the  major  region  of  the  surface  contributing  to  the  received  energy.  All 
values  below  -90  dB  have  been  set  equal  to  —90  dB. 
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Figure  14.  Scattering  strength  at  40  kHz  as  a  function  of  position  for  source  and 
receiver  positions  given  in  meters  at  left  (x,y,z),  Arafura  site.  There  is  no 
critical  angle  in  this  case.  The  angular  width  of  the  scattering-strength 
peak  level  (see  Figure  6)  is  so  narrow  that  higher  resolution  (Figure  15)  is 
needed  to  resolve  the  structure  above  the  receiver. 


Receiver 


Figure  15.  A  high-resolution  view  of  the  boxed  region  in  Figure  14. 
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Figure  1 8.  Scattering  strength  at  40  kHz  as  a  function  of  position  for  source  and 
receiver  positions  given  in  meters  at  left  (x,y,z),  Arafura  site. 
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Figure  19.  Scattering  strength  of  Figure  18  minus  (see  Eq.  7)  times  distance 
traveled  from  each  surface  position  to  the  receiver.  The  result  indicates 
the  major  region  of  the  surface  contributing  to  the  received  energy.  All 
values  below  -90  dB  have  been  set  equal  to  -90  dB. 
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Figure  20.  Scattering  strength  at  40  kHz  as  a  function  of  position  for  source  and 
receiver  positions  given  in  meters  at  left  (x,y,z},  Arafura  site. 
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Figure  21.  Scattering  strength  of  Figure  20  minus  (see  Eq.  7)  times  distance 
traveled  from  each  surface  position  to  the  receiver.  The  result  indicates 
the  major  region  of  the  surface  contributing  to  the  received  energy.  All 
values  below  -90  dB  have  been  set  equal  to  -90  dB. 
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When  the  source  is  20  and  10  m  above  the  surface  (Figures  9  and  11),  the  refracted  path 
in  the  sediment  is  short  enough  that  it  contributes  to  the  received  signal.  However,  when  the 
source  is  at  5  m  above  the  surface  (Figure  13),  the  contribution  via  the  refracted  path  has  un¬ 
dergone  so  much  attenuation  as  to  be  unimportant  at  the  receiver.  The  numerical  results  re¬ 
flected  in  Figure  13  indicate  that  the  primary  physics  contributing  to  the  received  intensity  is 
diffraction,  due  to  surface  roughness,  as  calculated  via  surface  perturbation  theory. 

Figures  14,  18,  and  20  show  the  scattered  strength  as  a  function  of  position  for  source 
heights  of  20, 10,  and  5  m,  respectively,  for  the  Arafura  site.  For  this  site,  the  elliptical  regions 
of  maximum  intensity  scattered  into  the  sediment  are  much  narrower  (see  Figures  5-7  and  the 
discussion  in  Section  4.1).  In  the  case  of  Figure  14,  in  fact,  a  higher-resolution  plot  (Figure  15) 
is  needed  to  see  the  region  in  the  vicinity  of  the  receiver.  Again,  it  is  useful  to  compare  the 
scattered  strengths  (Figures  14,  18,  and  20)  with  the  contributions  to  the  intensity  seen  at  the 
receiver  (Figures  16,  19,  and  21).  In  all  cases,  there  is  a  significant  contribution  along  the  re¬ 
fracted  path.  For  this  site,  the  sound  speed  in  the  sediment  is  less  than  that  in  the  water,  so  no 
critical  angle  exists.  As  in  the  Quinault  case,  much  of  the  received  intensity  comes  from  above 
the  receiver.  In  this  case,  however,  the  received  intensity  is  due  to  scattering  from  volume  in¬ 
homogeneities,  as  calculated  via  volume-perturbation  theory. 

4.3  Experimental  Tests  of  Scattering  Strength  Model 

Moe  et  al.  (1995)  used  the  surface-perturbation  model  of  Section  3.1.2  to  examine  the  ex¬ 
perimental  results  of  Chotiros  (1995).  The  results  could,  with  certain  assumptions  about  the 
surface  roughness  spectrum,  be  explained  using  the  model.  As  such,  Chotiros’s  experiment 
represents  a  consistency  test  of  at  least  one  aspect  of  the  present  model  with  real  data.  How¬ 
ever,  it  does  not  demonstrate  that  surface  perturbation  is  the  only  physics  contributing  to  the 
experimental  results.  Chotiros  has  shown  that  the  experimental  results  are  also  consistent  with 
a  Biot  slow  wave  propagating  in  the  sediment. 

The  geometries  of  Figures  8-21  were  chosen  with  an  experimental  test  of  this  model  in 
mind.  A  field  experiment  is  feasible  in  which  the  environmental  parameters  used  in  the  model 
and  the  intensity  at  buried  receivers  are  measured.  Obtaining  statistical  results  will  require  that 
the  source  in  the  water  be  movable.  An  experiment  of  this  type  has  recently  been  proposed. 
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5.  SUMMARY  AND  FUTURE  DIRECTIONS 

The  model  presented  here  for  scattering  into  the  sediment  includes  effects  of  both  surface 
roughness  and  volume  inhomogeneities.  It  closely  follows  the  development  of  Jackson  (1993). 
It  is  anticipated  that  in  the  long  term  the  Kirchhoff  and  surface-perturbation  theory  parts  of  the 
model  could  be  replaced  with  the  more  robust  small-slope  approximation  (Thorsos  and 
Broschat,  1995).  Likewise,  as  more  computational  power  becomes  available,  it  would  seem 
feasible  to  replace  the  treatment  of  the  volume  scattering  as  an  equivalent  surface  effect  with 
a  more  ambitious  calculation  of  scattering  within  the  sediment  volume. 

Having  said  this,  the  model,  as  developed  to  this  point,  seems  to  capture  many  of  the  ef¬ 
fects  one  might  anticipate  as  well  as  at  least  one  (subcritical  penetration)  that  might  not  be  so 
obvious.  Quantitative  testing  of  the  model  awaits  further  experiments. 
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APPENDIX 

Definition  of  Bistatic  Cross  Section 


Typically  the  bistatic  cross  section  is  defined  as 


a  = 


(Al) 


where  r  is  the  range  from  the  scattering  point  on  the  surface  to  the  receiver,  7^  is  the  scattered 
intensity  measured  at  r,  A  is  the  patch  of  the  surface  contributing  to  ,  and  /•  is  the  incident 
intensity.  However,  in  the  present  case,  the  effect  of  attenuation  in  the  sediment  is  taken  into 
account  in  addition  to  the  effect  of  spherical  spreading.  Thus  the  definition  of  cross  section 
used  here  is 


a  = 


^  2bk.r 


AL 


(A2) 


where  the  factor  in  parenthesis  accounts  for  the  attenuation  in  the  sediment,  is  the  real 
part  of  the  acoustic  wavenumber  in  the  sediment  (see  text  below  Eq.  3),  and  5  is  defined  in 
Table  1 .  This  definition  then  leads  to  correct  calculations  of  intensity  received  in  the  sediment 
when  attenuation  is  included  in  that  calculation. 
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